home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dposv.z / dposv
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPOOOOSSSSVVVV((((3333FFFF))))                                                            DDDDPPPPOOOOSSSSVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPOSV - compute the solution to a real system of linear equations  A * X
  10.      = B,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DPOSV( UPLO, N, NRHS, A, LDA, B, LDB, INFO )
  14.  
  15.          CHARACTER     UPLO
  16.  
  17.          INTEGER       INFO, LDA, LDB, N, NRHS
  18.  
  19.          DOUBLE        PRECISION A( LDA, * ), B( LDB, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      DPOSV computes the solution to a real system of linear equations
  23.         A * X = B, where A is an N-by-N symmetric positive definite matrix and
  24.      X and B are N-by-NRHS matrices.
  25.  
  26.      The Cholesky decomposition is used to factor A as
  27.         A = U**T* U,  if UPLO = 'U', or
  28.         A = L * L**T,  if UPLO = 'L',
  29.      where U is an upper triangular matrix and L is a lower triangular matrix.
  30.      The factored form of A is then used to solve the system of equations A *
  31.      X = B.
  32.  
  33.  
  34. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  35.      UPLO    (input) CHARACTER*1
  36.              = 'U':  Upper triangle of A is stored;
  37.              = 'L':  Lower triangle of A is stored.
  38.  
  39.      N       (input) INTEGER
  40.              The number of linear equations, i.e., the order of the matrix A.
  41.              N >= 0.
  42.  
  43.      NRHS    (input) INTEGER
  44.              The number of right hand sides, i.e., the number of columns of
  45.              the matrix B.  NRHS >= 0.
  46.  
  47.      A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
  48.              On entry, the symmetric matrix A.  If UPLO = 'U', the leading N-
  49.              by-N upper triangular part of A contains the upper triangular
  50.              part of the matrix A, and the strictly lower triangular part of A
  51.              is not referenced.  If UPLO = 'L', the leading N-by-N lower
  52.              triangular part of A contains the lower triangular part of the
  53.              matrix A, and the strictly upper triangular part of A is not
  54.              referenced.
  55.  
  56.              On exit, if INFO = 0, the factor U or L from the Cholesky
  57.              factorization A = U**T*U or A = L*L**T.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPOOOOSSSSVVVV((((3333FFFF))))                                                            DDDDPPPPOOOOSSSSVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDA     (input) INTEGER
  75.              The leading dimension of the array A.  LDA >= max(1,N).
  76.  
  77.      B       (input/output) DOUBLE PRECISION array, dimension (LDB,NRHS)
  78.              On entry, the N-by-NRHS right hand side matrix B.  On exit, if
  79.              INFO = 0, the N-by-NRHS solution matrix X.
  80.  
  81.      LDB     (input) INTEGER
  82.              The leading dimension of the array B.  LDB >= max(1,N).
  83.  
  84.      INFO    (output) INTEGER
  85.              = 0:  successful exit
  86.              < 0:  if INFO = -i, the i-th argument had an illegal value
  87.              > 0:  if INFO = i, the leading minor of order i of A is not
  88.              positive definite, so the factorization could not be completed,
  89.              and the solution has not been computed.
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.